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Abstract — In this paper, we study the performance of space 
modulation for Multiple-Input-Multiple-Output (MIMO) wire- 
less systems with imperfect channel knowledge at the receiver. 
We focus our attention on two transmission technologies, which 
are the building blocks of space modulation: i) Space Shift Key- 
ing (SSK) modulation; and ii) Time-Orthogonal-Signal-Design 
(TOSD-) SSK modulation, which is an improved version of SSK 
modulation providing transmit-diversity. We develop a single- 
integral closed-form analytical framework to compute the Aver- 
age Bit Error Probability (ABEP) of a mismatched detector for 
both SSK and TOSD-SSK modulations. The framework exploits 
the theory of quadratic-forms in conditional complex Gaussian 
Random Variables (RVs) along with the Gil-Pelaez inversion 
theorem. The analytical model is very general and can be used for 
arbitrary transmit- and receive-antennas, fading distributions, 
fading spatial correlations, and training pilots. The analytical 
derivation is substantiated through Monte Carlo simulations, and 
it is shown, over independent and identically distributed (i.i.d.) 
Rayleigh fading channels, that SSK modulation is as robust as 
single-antenna systems to imperfect channel knowledge, and that 
TOSD-SSK modulation is more robust to channel estimation 
errors than the Alamouti scheme. Furthermore, it is pointed 
out that only few training pilots are needed to get reliable 
enough channel estimates for data detection, and that transmit- 
and receive-diversity of SSK and TOSD-SSK modulations are 
preserved even with imperfect channel knowledge. 

Index Terms — Imperfect channel knowledge, "massive" 
multiple-input-multiple-output (MIMO) systems, mismatched 
receiver, performance analysis, single-RF MIMO design, space 
shift keying (SSK) modulation, spatial modulation (SM), 
transmit-diversity. 



I. Introduction 

SPACE modulation [1] is a novel digital modulation con- 
cept for Multiple-Input-Multiple-Output (MIMO) wire- 
less systems, which is receiving a growing attention due to 
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the possibility of realizing low-complexity and spectrally- 
efficient MIMO implementations [2]-[7]. The space modu- 
lation principle is known in the literature in various forms, 
such as Information-Guided Channel Hopping (IGCH) [2], 
Spatial Modulation (SM) [3], and Space Shift Keying (SSK) 
modulation [4] . Although different from one another, all these 
transmission technologies share the same fundamental working 
principle, which makes them unique with respect to conven- 
tional modulation schemes: they encode part of the informa- 
tion bits into the spatial positions of the transmit-antennas 
in the antenna-array, which plays the role of a constella- 
tion diagram (the so-called "spatial-constellation diagram") 
for data modulation [1], [7]. In particular, SSK modulation 
exploits only the spatial-constellation diagram for data mod- 
ulation, which results in a very low-complexity modulation 
concept for MIMO systems [4]. Recently, improved space 
modulation schemes that can achieve a transmit-diversity gain 
have been proposed in [8] — [1 1]. Furthermore, a unified MIMO 
architecture based on the SSK modulation principle has been 
introduced in [12]. 

In SSK modulation, blocks of information bits are mapped 
into the index of a single transmit-antenna, which is switched 
on for data transmission while all the other antennas radiate 
no power [4]. SSK modulation exploits the location-specific 
property of the wireless channel for data modulation [7]: 
the messages sent by the transmitter can be decoded at the 
destination since the receiver sees a different Channel Impulse 
Response (CIR) on any transmit-to-receive wireless link. 
In [4] and [7], it has been shown that the CIRs are the 
points of the spatial-constellation diagram, and that the Bit 
Error Probability (BEP) depends on the distance among these 
points. Recent results have shown that, if the receiver has 
Perfect Channel State Information (P-CSI), space modulation 
can provide better performance than conventional modulation 
schemes with similar complexity [2]-[4], [9], [10], and [13]- 
[17]. However, due to its inherent working principle, the major 
criticism about the adoption of SSK modulation in realistic 
propagation environments is its robustness to the imperfect 
knowledge of the wireless channel at the receiver. In particular, 
it is often argued that space modulation is more sensitive to 
channel estimation errors than conventional systems. The main 
contribution of this paper is to shed light on this matter. 

Some research works on the performance of space modu- 
lation with imperfect channel knowledge are available in the 
literature. However, they are insufficient and only based on 
numerical simulations. In [4], the authors have studied the 
ABEP of SSK modulation with non-ideal channel knowledge. 
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However, there are four limitations in this paper: i) the ABEP 
is obtained only through Monte Carlo simulations, which 
is not very much insightful; ii) the arguments in [4] are 
applicable only to Gaussian fading channels and do not take 
into account the cross-product between channel estimation 
error and Additive White Gaussian Noise (AWGN) at the 
receiver; iii) it is unclear from [4] how the ABEP changes 
with the pilot symbols used by the channel estimator; and 
iv) the robustness/weakness of SSK modulation with respect 
to conventional modulation schemes is not analyzed. In [18], 
we have studied the performance of SSK modulation when 
the receiver does not exploit for data detection the knowledge 
of the phase of the channel gains (semi-blind receiver). It 
is shown that semi-blind receivers are much worse than 
coherent detection schemes, and, thus, that the assessment of 
the performance of coherent detection with imperfect channel 
knowledge is a crucial aspect for SSK modulation. A very 
interesting study has been recently conducted in [19], where 
the authors have compared the performance of SM and V- 
BLAST (Vertical Bell Laboratories Layered Space-Time) [20] 
schemes with practical channel estimates. It is shown that the 
claimed sensitivity of space modulation to channel estimation 
errors is simply a misconception and that, on the contrary, 
SM is more robust than V-BLAST to imperfections on the 
channel estimates, and that less training is, in general, needed. 
However, the study in [19] is conducted only through Monte 
Carlo simulations, which do not give too much insights for 
performance analysis and system optimization. Finally, in [10] 
the authors have proposed a Differential Space-Time Shift 
Keying (DSTSK) scheme, which is based on the Cayley uni- 
tary transform theory. The DSTSK scheme requires no channel 
estimation at the receiver, but incurs in a 3dB performance 
loss with respect to coherent detection. Furthermore, it can be 
applied to only real-valued signal constellations. Unlike [10], 
which avoids channel estimation, we are interested in studying 
the training overhead that is needed for channel estimation 
and to achieve close-to-optimal performance with coherent 
detection. 

Motivated by these considerations, this paper is aimed at 
developing a very general analytical framework to assess 
the performance of space modulation with coherent detection 
and practical channel estimates. In particular, we focus our 
attention on two transmission technologies, which are the 
building blocks of space modulation: i) Space Shift Key- 
ing (SSK) modulation [4]; and ii) Time-Orthogonal-Signal- 
Design (TOSD-) SSK modulation, which is an improved ver- 
sion of SSK modulation providing transmit-diversity [8], [11]. 
Our theoretical and numerical results corroborate the findings 
in [19], and highlight three important outcomes: i) SSK 
modulation is as robust as single-antenna systems to imperfect 
channel knowledge; ii) TOSD-SSK modulation is more robust 
to channel estimation errors than the Alamouti scheme [21]; 
and iii) only few training pilots are needed to get reliable 
enough channel estimates for data detection. More precisely, 
we provide the following contributions: i) we develop a single- 
integral closed-form analytical framework to compute the 
Average BEP (ABEP) of a mismatched detector [22] for SSK 
and TOSD-SSK modulations, which can be used for arbitrary 



transmit- and receive-antennas, fading distributions, fading 
spatial correlations, and training pilots for channel estimation. 
It is shown that the mismatched detector of SSK and TOSD- 
SSK modulations can be cast in terms of a quadratic-form in 
complex Gaussian Random Variables (RVs) when conditioning 
upon fading channel statistics, and that the ABEP can be 
computed by exploiting the Gil-Pelaez inversion theorem [23]; 
ii) over independent and identically distributed (i.i.d.) Rayleigh 
fading channels, we show that SSK modulation is superior 
to Quadrature Amplitude Modulation (QAM), regardless of 
the number of training pulses, if the spectral efficiency is 
greater than 2 bpcu (bits per channel use) and the receiver 
has at least two antennas; iii) in the same fading channel, we 
show that TOSD-SSK modulation is superior, regardless of 
the number of antennas at the receiver and training pulses, to 
the Alamouti scheme with QAM if the spectral efficiency is 
greater than 2 bpcu. Also, unlike the P-CSI setup, TOSD- 
SSK modulation can outperform the Alamouti scheme if the 
spectral efficiency is 2 bpcu, just one pilot pulse for channel 
estimation is used, and the detector is equipped with at least 
two antennas; iv) still over i.i.d. Rayleigh fading, we show 
that, compared to the P-CSI scenario, SSK and TOSD-SSK 
modulations have a Signal-to-Noise-Ratio (SNR) penalty of 
approximately 3dB and 2dB when only one pilot pulse can be 
used for channel estimation, respectively. Also, single-antenna 
and Alamouti schemes have a SNR penalty of approximately 
3dB for QAM; and v) we verify that transmit- and receive- 
diversity of SSK and TOSD-SSK modulations are preserved 
even for a mismatched detector. 

The remainder of this paper is organized as follows. In 
Section [Til the system model is introduced. In Section [111] and 
Section |lVl SSK and TOSD-SSK modulations are described 
and the analytical frameworks to compute the ABEP with 
imperfect channel knowledge are developed, respectively. In 
Section [V] the spectral efficiency of TOSD-SSK modulation 
with time-orthogonal shaping filters is studied. In Section 
IVII numerical results are shown to substantiate the analytical 
derivation, and to compare SSK and TOSD-SSK modulations 
with state-of-the-art single-antenna and Alamouti schemes. 
Finally, Section IVHI concludes this paper. 

II. System Model 

We consider a generic N t x N r MIMO system, with N t 
and N r being the number of antennas at the transmitter and at 
the receiver, respectively. SSK and TOSD-SSK modulations 
work as follows [4], [8], [11]: i) the transmitter encodes blocks 
of log 2 (N t ) data bits into the index of a single transmit- 
antenna, which is switched on for data transmission while all 
the other antennas are kept silent; and ii) the receiver solves 
an A t -hypothesis detection problem to estimate the transmit- 
antenna that is not idle, which results in the estimation of the 
unique sequence of bits emitted by the encoder. With respect to 
SSK modulation [4], in TOSD-SSK modulation [11] the <-th 
transmit-antenna, when active, radiates a distinct pulse wave- 
form Wt(-) for t = 1, 2, . . . , Nt, and the waveforms across the 
antennas are time-orthogonal, i.e.0, J^°° Wt 1 (£) w* 2 (£) d£ = 

'(•)* denotes complex-conjugate. 
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if ti ± h and w tl (£) < (0 d£ = 1 if ti = h. In 
other words, SSK modulation is a special case of TOSD-SSK 
modulation with w t (£) = wq (£) for t = 1,2, ... , N t . 

In [11], we have analytically proved that the diversity order 
of SSK modulation is N r , while the diversity order of TOSD- 
SSK modulation is 2N r , which results in a transmit-diversity 
equal to 2 and a receive-diversity equal to N r . Thus, TOSD- 
SSK modulation provides a full-diversity-achieving (i.e., the 
diversity gain is N t N r ) system if Nt = 2. This scheme 
has been recently generalized in [24] to achieve arbitrary 
transmit-diversity. It is worth emphasizing that in TOSD-SSK 
modulation a single-antenna is active for data transmission 
and that the information bits are still encoded into the index 
of the transmit-antenna, and are not encoded into the impulse 
(time) response, w t (•), of the shaping filter. In other words, 
TOSD-SSK modulation is different from conventional Single- 
Input-Single-Output (SISO) schemes with Orthogonal Pulse 
Shape Modulation (O-PSM) [25], which are unable to achieve 
transmit-diversity as only a single wireless link is exploited 
for communication [11]. Also, TOSD-SSK modulation is 
different from conventional transmit-diversity schemes [26], 
and requires no extra time-slots for transmit-diversity. Further 
details are available in [11] and are here omitted to avoid 
repetitions. 

Throughout this paper, the block of information bits en- 
coded into the index of the i-th transmit-antenna is called 
"message", and it is denoted by m t for £ = 1,2,..., Nt. The 
Nt messages are assumed to be equiprobable. Moreover, the 
related transmitted signal is denoted by St (•). It is implicitly 
assumed with this notation that, if m t is transmitted, the analog 
signal Sf (•) is emitted by the t-th transmit-antenna while the 
other antennas radiate no power. 

A. Notation 

Main notation is as follows, i) We adopt a complex- 
envelope signal representation, ii) j = y— 1 is the imagi- 
nary unit, iii) (x <£> y) (u) — f_°° x (£) y (u ~ £) d£ is the 
convolution of signals x (•) and y(-). iv) |-| 2 is the square 
absolute value, v) E{-} is the expectation operator computed 
over channel fading statistics, vi) Re{-} and Im{-} are real 
and imaginary part operators, respectively, vii) Pr{-} denotes 
probability, viii) Q (u) = (l/V^F) J u + °° exp (-£ 2 / 2 ) d(, is 
the Q-function. ix) 8 (•) and <5. . are Dirac and Kronecker 
delta functions, respectively, x) Mx (s) — E {exp (sX)} and 
{v) — E {exp (juX)} are Moment Generating Function 
(MGF) and Characteristic Function (CF) of RV X, respec- 
tively, xi) cx denotes "is proportional to". 

B. Channel Model 

We consider a general frequency-flat slowly-varying chan- 
nel model with generically correlated and non-identically 
distributed fading gains. In particular (t — 1,2,..., Nt, r = 

1,2 Y,..,: 

• ht,r (0 = a t.r5 (£ — Tt ir ) is the channel impulse re- 
sponse of the transmit-to-receive wireless link from 
the t-th transmit-antenna to the r-th receive-antenna. 



ot t ,r = Pt,r ex p(J l Pt,r) i s me complex channel gain with 
Pt r and (p t r denoting the channel envelope and phase, 
respectively, and r t r is the propagation time-delay. 
• The time-delays Tt, r are assumed to be known at the 
receiver, i.e., perfect time-synchronization is considered. 
Furthermore, we consider ri,i = T1.2 — ■ • • = tn u n t , 
which is a realistic assumption when the distance between 
the transmitter and the receiver is much larger than the 
spacing between the transmit- and receive-antennas [7]. 
Due to the these assumptions, the propagation delays can 
be neglected in the remainder of this paper. 

C. Channel Estimation 

Let E p and N p be the energy transmitted for each pilot pulse 
and the number of pilot pulses used for channel estimation, 
respectively. Similar to [27] and [28], we assume that channel 
estimation is performed by using a Maximum-Likelihood 
(ML) detector, and by observing N p pilot pulses that are 
transmitted before the modulated data. During the transmission 
of one block of pilot-plus-data symbols, the wireless channel 
is assumed to be constant, i.e. a quasi-static channel model is 
considered. With these assumptions, the estimates of channel 
gains a.tr (t = 1,2,..., N t , r = 1,2,..., N r ) can be written 
as follows: 

a t ,r = Pt,r exp (J0t,r) = Oi t ,r + St,r (1) 

where &t,r, $t.r, and ipt, r are the estimates of f3 t ,r, and 
<fit,r> respectively, at the output of the channel estimation unit, 
and e t r is the additive channel estimation error, which can be 
shown to be complex Gaussian distributed with zero-mean and 
variance = N /(E p N p ) per dimension [27], [28], where 
iVo denotes the power spectral density per dimension of the 
AWGN at the receiver. The channel estimation errors, e t r , are 
statistically independent and identically distributed, as well as 
statistically independent of the channel gains and the AWGN 
at the receiver. 

D. Mismatched ML-Optimum Detector 

For data detection, we consider the so-called mismatched 
ML-optimum receiver according to the definition given in 
[22]. In particular, a detector with mismatched metric estimates 
the complex channel gains as in (Q~|), and uses them in the same 
metric that would be applied if the channels were perfectly 
known. To avoid repetitions in the analysis of SSK and TOSD- 
SSK modulations, the mismatched detector is here described 
by assuming arbitrary shaping filters. 

The mismatched ML-optimum detector can be obtained as 
follows. Let m q with q = l,2,...,N t be the transmitted 
message. The signal received after propagation through the 
wireless fading channel and impinging upon the r-th receive- 
antenna can be written as follows: 

7 -r (0 = s qyr (£) + r} r (£) if m q is sent (2) 

where: i) s 9 , r (£) = (s q (g) h q , r ) (£) = a q ^ r s q (£) = 
Pq, r exp(j<Pq,r)sq(ti) for q = 1, 2, . . . , N t and r = 
1,2,...,AV; ii) s q (£) = VE^w q (0 for q = 1, 2, . . . , N t , 
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D mq (m t ) = -Y] \ 

r=l LJt 



E 



(4) 



-Dm, (m t ) 



E 

r=l 



ira?,rMo (4) 



Vr(0 



m i I Em , & \ 

— a t} rW (£) + J —St,rWQ (£) 

iVo V -^0 



d£ > (5) 



?7i = arg max 

m t for t=l,2,...,JV, 



{-Dm, (m t )} 



oc arg mm 

m t for t=l,2,...,iV t 



{^K (m*)} 



arg mm 

m t for t=l,2,...,JV, 



VO,; 



— (a tir -a qir ) + ] j—e t>r 



(6) 



where i? m is the average energy transmitted by each an- 
tenna that emits a non-zero signal; and iii) rj r (•) is the 
complex AWGN at the input of the r-th receive-antenna for 
r = 1,2, . . . , N r , which has power spectral density Nq per 
dimension. Across the receive-antennas, the noises r\ r (■) are 
statistically independent. 

In particular, (0 is a general iV t -hypothesis detection prob- 
lem [29, Sec. 7.1], [30, Sec. 4.2, pp. 257] in AWGN, when 
conditioning upon fading channel statistics. Accordingly, the 
mismatched ML-optimum detector with imperfect CSI at the 
receiver is as follows: 



m 



arg max 

mt for t=l,2,...,iVt 



{-Dm, (mt)} 



(3) 



where rh is the estimated message and D m (mt) is the mis- 
matched decision metric [7], [11], which is shown in (0 on top 
of this page, where s t . r (£) = ct t , r s t (£) = (a t , r + e t . r ) s t (£) 
and T m is the symbol period. 

III. SSK Modulation 

A. Decision Metrics 

In SSK modulation, the decision metric in (0 can be re- 
written from (0 and (0 as shown in on top of this page, 
where we have taken into account that for SSK modulation the 
shaping filters are all equal to wq (■), and we have introduced 
the scaling factor 1/Nq, which does not affect 0. 

From (0, and after some algebra, the maximization problem 
in (0 reduces to (0 shown on top of this page, where 
f)o,r = J Tm Vr (0 ™o (£) d£, and (m t ) is statistically 
equivalent to D m (mt). In particular, (0 can be thought as a 
mismatched detector in which: i) first, pulse-matched filtering 
is performed; and ii) then, ML-optimum decoding is applied 
to the resulting signal. 

B. ABEP 

The ABEP of the detector in (0 can be computed in closed- 
form as follows: 



n (mm . 

ABEp'i'E 

U=i t=i 1 



N H (t,q) 



log 2 (N t ) 



Pr { m — m t \ m q } 



(7) 



-j^JV t log a (N t )- 



(a) (6) 

where = comes from [31, Eq. (4) and Eq. (5)], and = 
is the asymptotically-tight union-bound recently introduced 
in [7, Eq. (35)]. Furthermore, Nu(t,q) is the Hamming 
distance between the bit-to-antenna-index mappings of mt 
and m q ; and APEP (m q — > m t ) = E {PEP (m q — > mt)} = 
E{Pr{m 9 — > m t }} is the Average Pairwise Error Probability 
(APEP), i.e., the probability of estimating m t when, instead, 
m q is transmitted, under the assumption that m t and m q are 
the only two messages possibly being transmitted. 

Let us note that (0 simplifies significantly when 
APEP (TO, -> mt) = APEP for t = 1,2,..., N t and 
q = 1, 2, . . . , Nt (e.g., for i.i.d. fading). In this case, the ABEP 
in becomes: 

ABEP < J^r j:j:N H (t,q)^^APEP (8) 



NtlogAN)^ 



(a) 



where = comes from the identity E 9 =i St=i Nh (t, q) = 
(Nt /2) log 2 (Nt), which can be derived via direct inspection 
for all possible bit-to-antenna-index mappings. 



C. Computation of PEPs 

Let us start by computing the PEPs, i.e., the pairwise 
probabilities in (0 when conditioning upon fading channel 
statistics. From (0, PEP (m q — > mt) is as follows: 



PEP 
where: 

,(e) 



(m q -> m t ) = Pr (m t ) < (m,)} (9) 



D)n' q (m t ) = E 
r=l 

^! (m,) = E 



>?0,r 



10,r 

■Sm 




By introducing the notation (r = 1,2,..., N r ): 




vm 

VO.r 

vm 



(at; 



Nq £ 1, r 



(ID 



the PEP in (0 can be re-written in the general form (with 

A = 1, B = -1, and C = 0): 



APEP(m,->mi) 



PEP (m q -> mt) = Pr {D < 0} 



(12) 
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*D (H at,q) = - N — ]v7 



exp 



(-V 2 galt,q+ jV9blt,q) 

(y + jv a ) (y - jv b ) 



THcxp{A(^ 7M } (15) 



ABEP < 




(19) 



where: 



£> = ^ U |A r | 2 + B \Y r \ 2 + CX r r; + C*X*,Y r j (13) 

r=l 

From [23, Sec. Ill], we notice that, when conditioning upon 
fading channel statistics, the RV D in ( fT3l ) is a quadratic-form 
in complex Gaussian RVs. In fact, AWGN at the receiver input 
and channel estimation error are Gaussian distributed RVs. 
Furthermore, they are mutually independent among themselves 
and across the N r receive-antennas. Literature on quadratic- 
forms in complex Gaussian RVs is very rich, and during the 
last decades many different techniques have been developed 
for their analysis (see, e.g., [23] and [32] for a survey). 
Furthermore, effective methods for the computation of the PEP 
over generalized fading channels have been proposed, e.g., 
[33]-[35], and simple analytical frameworks for some special 
fading scenarios are available in [36, Ch. 9]. In this paper, we 
propose to use the Gil-Pelaez inversion theorem [37]. 

Accordingly, by using [37] the PEP can be computed as 
follows: 



PEP (m, 



m t ) 



1 1 

2 ~n jQ 

1 1 

2 ~ ~ 



+ 00 



Im D (v\ a t , ,)} 



dv 



7T 



^ 2 Im{tt g (tan(fl|a t , q )} 



o 



sin (0 cos (£) 



(14) 



where typ ( • | ott, q ) is the CF of RV D when conditioning upon 



the channel gains, and a. t , q = {a t , r , a q,r} r li I s a short-hand 
to denote all the channel gains in ( TT3T >. 

The conditional CF, ^ d ('\ ctt.q), of RV D is given by 
[23, Eq. (2) and Eq. (3)] shown in ( fTBT l on top of this 
page, where: i) j=E m /N ; ii) r pm = E p /E m ; iii) g a = 
27 1 + {Nprpmy 1 ; iv) g b = 7; and: 



r„ = v b = (1/2) \l (N p r pm ) 2 + 2{N p r pm ) 1 



N,. 



Jt,q =7 («*,?) = E 



r=l 



at,: 



A.(i/) = v a v b {-v 2 g a + jvg b ) (v + jv a ) (y - jv b ) 1 

(16) 



, T(v) = (v a v b ) Nr (v + jv a ) Nr (v - jvb) Nr 



D. Computation of APEPs 

The APEP can be computed from (fT4l by removing the 
conditioning over the fading channel: 



APEP (m q -» m t ) = E {PEP (r 

_ 1 

~ 2 ~ ' 



m t )} 



+ °° Imft D {v)) dv (17) 



where V D (v) = E{* D «*,«)} is the CF of RV D 
averaged over all fading channel statistics. It can be computed 
from (Q3), as follows: 

* D (y) = E {T (1/) exp {A ( V ) 7t ,J} ( = } T (1/) M 7t , 9 (A (y)) 

(18) 

, s (a) 

where M lt (•) is the MGF of RV 7 t i9 , and = comes from 
the definition of MGF. 

In conclusion, the ABEP of SSK modulation over arbitrary 
fading channels and with practical channel estimates can be 
computed in closed-form from (0, (fTTI i, and ( fT8l , as shown 
in ( fl9] l on top of this page. 

The formula in (|T~9T > provides a very simple analytical tool 
for performance assessment of SSK modulation with channel 
estimation errors, and allows us to estimate the number of pilot 
pulses, N p , and the fraction of energy, r pm , to be allocated to 
each pilot pulse to get the desired performance. In particular, 
( fT9l needs only the MGF of RV 7t g to be computed. This 
latter MGF is the building block for computing the ABEP 
with P-CSI, and it has been recently computed in closed- 
form for a number of MIMO setups and fading conditions. 
In particular: i) it is known in closed-form for arbitrary 
correlated Nakagami-m fading channels and N r = 1 [7]; ii) it 
can be derived from [7] for independent Nakagami-m fading 
channels and arbitrary N r [38]; iii) it can be derived from 
[7] for Nakagami-m fading channels and arbitrary N r when 
the channel gains are correlated at the transmitter-side but 
are independent at the receiver-side [38]; and iv) it is known 
in closed-form for arbitrary correlated Rician fading channels 
and arbitrary N r [11]. For example, for i.i.d. Rayleigh fading 
channels, the ABEP in ( fl9] > reduces, from (0, to: 



ABEP < — - — 

~ 4 



2n Jo 



1 



(i-2n„A(i/)) 



'(20) 



where Slo = E ||a t J . | 2 | is the mean square value of the i.i.d. 
channel gains. 

Finally, we conclude this section with three general com- 
ments about ( fT9l : i) the integrand function is, in general, well- 
behaved when v ~ > for typical MGFs used in wireless 
communication problems. Thus, the numerical computation of 
the integral does not provide any critical issues. The interested 
reader might check this out in ( |20l , where it can be shown that 
the integrand function tends to a finite value when v — > 0; ii) 
since the ABEP depends on the MGF of RV j tt q, from [11] 
and [39] we conclude that the diversity order of the system is 
given by N r , which is the same as the P-CSI scenario. We 
will verify this statement in Section [VTl with some numerical 
examples, which will highlight that there is no loss in the 
diversity order with practical channel estimation; and iii) by 
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N r N r 

Dm q (mt) = Re {{ a t,r + £t,r)* VE^iVt.r} ~ E \ & t,r\ 

r—1 r—1 

N r N r 

Dm, K) = R-e{(a 9 , r + £g,r)* {a q , r E m + s/E^fjq^)} ^ I <*«,»• I 

r—1 r—1 



(22) 



PEP (m, -> m t ) 

N r 



- Pr < 



E [l ("9,rv / 7 + £ 9,''V / 7)* ( Q 9' r V / 7 + + I ( a <7,r V / 7 + £ 9,'-V / 7) (<*«,r- V / 7 + ~ | K,r + £ 9 ,r ^7 f] 



(24) 



E 

r=l 



direct inspection, it can be shown that dl9l ) reduces to the 
P-CSI lower-bound if N p r pm —> +00. 

IV. TO SD-SSK Modulation 

In this section, we focus our attention only on decision 
metrics and PEPs/APEPs since (Q and ([8]) are general and 
can be used for TOSD-SSK modulation too. 

A. Decision Metrics 

In TOSD-SSK modulation, the decision metric in ©, can 
be re-written as: 

N r 



D mq {mi) = ^ Re ^a q . r a* t r E m 8 t . q 



E m ^ — • , 2 



r=l 

(21) 

where we have taken into account that the shaping filters, 
w t (•), are time-orthogonal to one another, and we have 
defined f\ t r = J T rj r (£) w* t (£) <i£. In particular, for t ^ 9 
and t — q, the decision metric in (1211 ) simplifies as shown in 
4221 on top of this page. 

B. Computation of PEP s 

The PEPs, PEP (m q — > mi), in (O can be computed from 
© and 423: 



PEP (m 9 m t ) = Pr { D m? (m ? ) < £> TOg (m t )| 
= Pr 



D mq (m q ) D„ lq {mi) \ ( 23 ) 



No 



< 



N 



By using the identity Re{ab*} = (1/2) ab* + (1/2) 0*6, 
which holds for every pair of complex numbers a and 6, and 
by explicitly showing the SNR *f=E m /N in 423, the PEPs 
in 4231 simplifies as shown in d24i l on top of this page. 
By introducing the RVs (r = 1, 2, . . . , N r ): i) X q ^ r = 



0>-q,r\fl + £ q ,r\fl\ ») Yq, 



q , r 



V7 + {Vq,r/VNa); iii) 



■Xf.r = at,ry/l + et,rV7; iv) F t , r = rjt,r / V^O", an & 

D q =T. (A\X q , r \ 2 +B\Y q , r \ 2 + CX q ,rY*+C*X* Y q , 
r=l V 

a= E (A|Xi, r | 2 + B|yt, r | 2 + c7Xt, r i;* r + cfx t * r yt, r 

r=l v 



(26) 



the PEP in 42H» can be re-written (with A = -1/2, B = 0, 
and C = 1/2) as: 

PEP {m q -> mt) = Pr{£>, < A} 

= Pr{A, 9 = D, - A < 0} 

Similar to Section UlI-Cl from [23, Sec. Ill], we can readily 
conclude that both D q and D t in 425l l are quadratic-forms 
in conditional complex Gaussian RVs. Furthermore, we note 
that D q and D t are, when conditioning upon the fading 
channel gains, statistically independent, as AWGN and channel 
estimation errors are independent from one another if t 7^ q 
and across the N r receive-antennas. We emphasize that to 
compute the ABEP in (O we are interested only in the cases 
where t 7^ q, as Nh {t, q) = if t = q. 

From 4261 ). the PEPs can be still computed by using the 
Gil-Pelaez inversion theorem [37]: 

1 1 f+°° lm{*D t , q {v\a t , q )} 
2 



PEP {m q -> m t ) 



-dv 



n Jo v 
(0)1 _ 1 f+°° Im{V Dq (v\a q )V Dt (-v\att)} 
2 : 



dv 
(27) 



where ^ jj t q {-\ ot t q ) is the CF of RV D t q when con- 
ditioning upon the fading gains ot t , q = {(*t,r, a g,r} r= Ii, 
and *d„ ('|ag) and * Dt (-|a t ) are the CFs of RVs D q 
and D t when conditioning upon the fading gains a q = 
{a?,r}^ 1 and a t = {a*,r} r= Ii, respectively. Furthermore, 

(a) 

= comes from the independence of the conditional RVs D q 
and Dt, and the definition of CF, i.e., ^D t . (v\ ctt, q ) = 
Ej7,e {exp (jvDt, q )} = E^, £ {exp(jVDg)exp(-.jVD i )} = 
i&D {v\ Q-q) ^D t {— v\ at-i). We emphasize that Fi v , e {■} is 
the expectation operator computed over AWGN and channel 
estimation errors, as we are conditioning upon the channel 
gains. 

The last step is to compute the CFs in d27l >. which can be 
obtained from [23, Eq. (2) and Eq.(3)] by using the theory of 
quadratic-forms in conditional complex Gaussian RVs, as: 

/ *D,(^|a 9 ) =T 9 (i/)exp{A 9 (j/)7 g } 

\ * Dt act) = T t (1/) exp { A t (1/) 7t } U8j 



where we have defined: i) u a 

ii) v b 



VW4) 



(1/2); iii) 7, 



|a 9 ,r| 2 ; iv) 7t = 7(a { ) = E^i l«*^| 2 ; v) ff a ' 



^(1/4) + N p r prn + (1/2); 

7(«g) = 
(9) _ 



(25) 



(1/2)7 l + ^v,,.)- 1 ; vi) - «M - 



9b =9a 



■^ = (1/2)7; 
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ABEP < 



N t log 2 (N t ) 



EE 

q=l t=l 



v»«/.,)|i-i 



+00 Im 



dv 



(32) 



ABEP<^-^r° InJ T « MT * ( -^ 
~ 4 2vr 



f (1 - n Q A q (v)f r (1 - ft A t (-*/)) 



(34) 



and: 



A, (*) = „ „ 6 + jvg[ q) ) (y + jv^- 1 {v - j^)" 1 

A t (f) = v a v b (-v 2 g^ + jvg^J (v + jva)' 1 (u - jvt,)- 1 

r q („) = T t H = (e.^)"- (1/ + jv a )~ Nr (y - m)~ Nr 

(29) 

C. Computation of APEPs 

The APEP can be computed from (|27| and ([28j by still 
using the Gil-Pelaez inversion theorem [37]: 

APEP (m q -> m t ) = E {PEP (m, -> to*)} 



+°°Im{vI/ Dt »} ^ (30) 
-dv 



where ®D t , q {v) = E{^ Dg (v\a q )^n t {-v\ot t )}, which, 
for generic fading channels, is: 

*J>*.« W = E 01 a,) *A (-H a*)} 

= E {T q (u) exp { A q (1/) 7 J T t (-1/) exp { A t (-1/) 7 t}} 
= T g (1/) T 4 (-1/) E {exp {A g {v) lq + A t (-u) lt }} 
= T q {v) T t (-i/)M 7 (A, w (l) 

(31) 



E 



is the MGF of RV 



and M (A), . (s) 

(i/)=A g (i/) 7 , + A t (-i/) 7t . 

In conclusion, the ABEP of TOSD-SSK modulation over 
arbitrary fading channels and with practical channel estimates 
can be computed in closed-form from (jTj, (l30l , and (l3TT l as 
shown in (f32l > on top of this page. 

Similar to SSK modulation, d32i > is general and useful for 
every MIMO setups. To be computed, a closed-form expres- 
sion of the MGF of RV 7t ( ^ } (v) = A q (v) j q + A t {-v) 7t , 
which is given by the linear combination of the power-sum of 
genetically correlated and distributed channel gains, is needed. 
This MGF is available for various fading channel models in 
[29], or, e.g., it can be readily computed by exploiting the 
Moschopoulos method for arbitrarily correlated and distributed 
Rician fading channels, as described in [11]. In particular, 
if all the channel gains are independent, but not necessarily 
identically distributed, the MGF M^ja)^ (•) reduces to: 

M 7 (A )M ( S ) = E{exp( S7t ( ^ (1/))} 

= E {exp {sA q (v) lq }} E {exp {sA t {-v) 7f }} 



l[M laqrl 2( s A q („)) 



r=l 



JJ M lat r]2 (sA t (-v)) 



r=l 



(33) 

where the MGFs M^ at ^2 (•) and Mj a r ^ (•) are available 



in closed-form in [29] for almost all fading channel models 
of interest in wireless communications. For example, if the 
channel gains are i.i.d. Rayleigh distributed the ABEP in (|32] > 
reduces, from dD, to (l3~4l shown on top of this page. 

Finally, similar to Section IIII-DI we note that: i) the 
integrand function in (f32l > is, for typical MGFs used in 
communication problems, well-behaved when v — > 0; ii) since 
the ABEP in, e.g., (l33l and (l3~4l is given by the products of 
2N r MGFs, we conclude from [11] and [39] that the diversity 
order of the system is 2N r , which is the same as the P-CSI 
scenario [11]; and iii) (l32l reduces to the P-CSI lower-bound 
if N p r prn -> +00. 

V. Bandwidth Efficiency of Orthogonal Shaping 
Filters Design 

In Section lIVl we have shown that TOSD-SSK modulation 
provides, even in the presence of channel estimation errors and 
with a single active antenna at the transmitter, a diversity order 
that is equal to 2N r . This is achieved by using time-orthogonal 
shaping filters at the transmitter, which is an additional design 
constraint that might not be required by SSK modulation and 
conventional single- and multiple-antenna systems. Thus, for 
a fair comparison among the various modulation schemes, 
it is important to assess whether the time-orthogonal con- 
straint affects the overall bandwidth efficiency of the com- 
munication system. More specifically, this section is aimed 
at understanding whether a larger transmission bandwidth is 
required for the transmission of the same number of bits in a 
given signaling time-interval T m , i.e., for a given bit/symbol 
or bpcu requirement. To shed light on this matter, in this 
section we analyze the bandwidth occupancy of commonly 
used shaping filters and, as an illustrative example, a family 
of recently proposed spectrally-efficient orthogonal shaping 
filters. More specifically: i) as far state-of-the-art shaping 
filters are concerned, we consider well-known time-limited 
rectangular, half-sine, and raised-cosine prototypes [40, Sec. 
III-B]; on the other hand, ii) as far as time-orthogonal 
shaping filters are concerned, we consider waveforms built 
upon linear combinations of Hermite polynomials [11], [25]. 
The analytical expressions of time and frequency responses of 
these letter filters are available in Appendix U for Nt — 4. 

Three important comments are worth being made about the 
shaping filters that are considered in our comparative study: 
1) we limit our study to considering time-limited shaping 
filters, which are simpler to be implemented than bandwidth- 
limited filters [40], [41]. This choice allows us to perform 
a fair comparison among SSK modulation and conventional 
modulation schemes. In fact, an important benefit of SSK and 
TOSD-SSK modulations is to take advantage of multiple- 
antenna technology with a single Radio Frequency (RF) front 
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TABLE I 

Bandwidth of various time-limited shaping filters . Time and frequency responses of rectangular, half-sine, and raised-cosine 

shaping filters are available in the captions of flg.fjj the shaping filters w t (•) are given in l |35l . let 
P (lu) = (l/y/^j f^^P (?) ex P ( — i^O d£, BE the Fourier transform of a generic shaping filter with time responsep (■). Then: i) the 

Fractional Power Containment Bandwidth (FPCB) is defined as FPCB X „, = min < B\ fe^r^-^-j- > X% I [48, p. 15], which is 

/o B6[0, + oo) I. So \P(u)\ 2 duj J 
THE BANDWIDTH B WHERE X % PERCENT OF THE ENERGY IS CONTAINED; AND II) THE BOUNDED POWER SPECTRAL DENSITY BANDWIDTH (BPSDB) 

IS DEFINED AS BPSDB TH(1B = min { B\ log 10 (\P(ui)\ 2 ) < log 10 (\P (^ pe ak)| 2 ) - TH dB , > B ) [48, P. 18], WHICH IS THE 
BANDWIDTH B BEYOND WHICH THE SPECTRAL DENSITY IS TH dB BELOW ITS PEAK (MAXIMUM VALUE), i.e., \P (w poak )| 2 . 



Fractional Power Containment Bandwidth (B/(2tt) kHz) 


x % 


Rectangular 


Half-Sine 


Raised-Cosine 


«*(•) 


99% 


7.61 


1.18 


1.41 


4.97 


99.995% 


>30 


6.98 


3.29 


6.46 


99.9999% 


>30 


22.14 


6.64 


7.31 


99.99999% 


>30 


29.96 


10.57 


7.76 


Bounded Power Spectral Density Bandwidth (B/(2w) kHz) 


TH dB 


Rectangular 


Half-Sine 


Raised-Cosine 


wt(-) 


3dB 


9.59 


2.28 


1.85 


6.35 


5dB 


>30 


8.18 


4.62 


7.40 


6dB 


>30 


15.13 


6.64 


7.85 


7dB 


>30 


28.03 


9.65 


8.27 


lOdB 


>30 


>30 


>30 


9.39 



end at the transmitter [3], [4], which is a research challenge 
that is currently stimulating the development of novel MIMO 
concepts based, e.g., on parasitic antenna architectures [42]- 
[44] . A recent survey on single-RF MIMO design is available 
in [45]. In order to use a single-RF chain, SSK and TOSD- 
SSK modulations need shaping filters that are time-limited 
and have a duration that is equal to the signaling time- 
interval T m . In fact, as remarked in [4, Section II-D], the 
adoption of shaping filters that are not time-limited would 
require a number of RF chains that is equal to the number 
of signaling time-intervals T m where the filter has a non- 
zero time response {i.e., the time-duration of the filter). Thus, 
bandwidth-limited shaping filters [41] would require multiple 
RF chains; 2) even though the orthogonal shaping filters 
considered in the present paper and summarized in Appendix 
Hare obtained by using the algorithm proposed in [25], which 
was introduced for Ultra Wide Band (UWB) systems, time- 
duration and bandwidth can be adequately scaled for narrow- 
band communication systems. For example, Fig.[T]is represen- 
tative of a narrow-band system with pulses having a practical 
time-duration of milliseconds and a practical bandwidth of 
kilohertz. Thus, neither UWB nor Spread Spectrum (SS) 
systems with orthogonal spreading codes are needed for space 
modulation; and 3) the method proposed in [25] for the design 
of orthogonal shaping filters guarantees that all the waveforms 
have the same time-duration and (practical) bandwidth. Thus, 
unlike conventional Hermite polynomials, time-orthogonality 
is guaranteed without bandwidth expansion. Let us emphasize 
that other methods are available in the literature to generate 



time-limited and time-orthogonal shaping filters. Two ex- 
amples, which allow us to jointly tuning time-duration and 
bandwidth and to guaranteeing low out-of-band interference, 
are given in [46] and [47]. 

Let us now compare the bandwidth efficiency of the orthog- 
onal shaping filters available in Appendix [Q with state-of-the- 
art shaping filters. A qualitative and quantitative comparisons 
are shown in Fig. [T] and in Table |U respectively, by using 
two commonly adopted definitions of bandwidth [48]: i) the 
Fractional Power Containment Bandwidth (FPCB) [48, p. 
15]; and ii) the Bounded Power Spectral Density Bandwidth 
(BPSDB) [48, p. 18]. The formal definition of these two 
concepts of bandwidth is given in the caption of Table [I] 
By carefully analyzing both Table U and Fig. [T] we notice 
that the bandwidth efficiency of the different shaping filters 
depend on how stringent the criterion to define the bandwidth 
is. In particular, if the percentage of energy that is required to 
be contained in the bandwidth (FPCB) is 99%, then the best 
shaping filter to use is the half-sine. On the other hand, if, 
to reduce the interference produced in adjacent transmission 
bands, the requirement moves from 99% to 99.99999%, then 
the best shaping filters to us are those given in Appendix U 
A similar comment applies when the BPSDB definition of 
bandwidth is used, but the best shaping filters are the raised- 
cosine (less stringent requirement) and the orthogonal filters in 
Appendix U (more stringent requirement). A similar trade-off 
has been shown in [40] and [48] for conventional modulation 
schemes and shaping filters. 

In other words, the shaping filters in Appendix U are 
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<o/(2n) [Hz] 



Fig. 1. Examples of time-limited shaping filters commonly used in 
the literature (frequency responses). The time responses are as follows, i) 
Rectangular pulse: v (£) = Pt (£)> where pt (£) = 1 if —To/2 < 
£ < To/2 and pt (£) = 1 elsewhere, ii) Half-sine pulse: v (£) = 
v2 sin [tt + 0.5To) /To] pt (£)■ iii) Raised-cosine pulse: u (f) = 
vWsjl - cos [2tt (£ + 0.5T ) /T ]}p To (£)• iv ) The orthogonal shaping 
filters for JV t = 4 are given in 051 in Appendix U with t = 1CT 4 . T = 
10 — 3 is the time duration of the filters. The frequency response is defined 
as V (uj) = (l/y/2wj f^^v(£)exp (—jw£)d£. They are as follows, i) 
Rectangular pulse: | V (2ttuj)\ = K r sinc (cjTq), where sine (a;) = 1 if X = 
and sine (a;) = sin (irx) /irx if x ^ 0, ii) Half-sine pulse: \V (27nj)| = 
re/js [cos (wajTo) / (l — 4w 2 Tq )] . iii) Raised-cosine pulse: |V(27T(d)| = 
Krc [sine (wT ) + (1/2) sine (uiTq - 1) + (1/2) sine (luTq + 1)]. iv) The 
orthogonal shaping filters for Nt = 4 are given in {36) in Appendix [I] with 
to = 10 -4 . K r , Kjis, and K rc are constant factors that are not relevant for 
our analysis. 

designed to have a very flat spectrum in the transmission band 
to improve the energy efficiency, as well as a very fast roll-off 
to reduce interference and enhance coexistence capabilities. 
This is especially useful to increase the system efficiency 
since current standards require the transmitted spectrum to 
occupy a well-defined spectral mask, e.g., for Wireless Local 
Area Networks (WLAN) and UWB wireless systems. For this 
reason, the shaping filters in Appendix U have a very good 
energy containment and bounded energy spectrum. Finally, we 
emphasize that the shaping filters given in Appendix Q] are just 
an example of time-orthogonal filters that can be obtained with 
state-of-the-art signal processing algorithms [46], [47], as 
well as that the waveforms compared in Table U have the same 
time-duration T m , and, thus, they provide the same signaling 
rate 1/T m . 

For illustrative purposes, in this paper we choose the 
shaping filters with the main objective to limit, as much 
as possible, out-of-band interference in order to enhance 
the coexistence capabilities of our communication system, 
and to reduce interference in adjacent transmission bands. 
Thus, our criterion is based on choosing filters which, for 
the same time-duration, have a stringent energy containment 
or bounded energy spectrum. For example, we assume either 
X% > 99.9999% or TH dB > 6dB in Table H With these 
assumptions, the orthogonal shaping filters given in Appendix 
U are the best choice, and are chosen to obtain the simulation 
results in Section |VI] For applications where less stringent 



coexistence capabilities might be required, the shaping filters 
given in Appendix [Q might not be the best choice, as they 
would require a larger bandwidth. In that case, by using the 
algorithms in [25], [46], [47], and references therein, new 
orthogonal pulses could be generated with the required time- 
duration and (practical) bandwidth. 

VI. Numerical and Simulation Results 

In this section, we show some numerical examples in order 
to: i) study the performance of SSK and TOSD-SSK modula- 
tions in the presence of channel estimation errors; ii) compare 
the achievable performance with single-antenna and Alamouti 
schemes; and iii) assess the accuracy of our analytical deriva- 
tion. For illustrative purposes, i.i.d Rayleigh fading channels 
are considered in all the analyzed scenarios. The interested 
reader might find in [7], [11], and [18] numerical examples 
about the performance of SSK and TOSD-SSK modulations 
for different wireless channels. Single-antenna and Alamouti 
schemes are chosen as state-of-the-art transmission technolo- 
gies for performance comparison because they have the same 
diversity order and the same decoding complexity as SSK and 
TOSD-SSK modulations, respectively. The interested reader 
might find in [49, Fig. 2] the comparison with transmit- 
diversity Space-Time-Block-Codes (STBCs) for MIMO sys- 
tems with more than two antennas at the transmitter, and in 
[50, Fig. 8] the comparison with spatial multiplexing MIMO 
systems with multi-user detection. In these latter cases, both 
STBCs and spatial multiplexing MIMO have higher decoding 
complexity and worse performance than space modulation. 

The simulation setup used in our study is as follows: i) we 
consider i.i.d Rayleigh fading with unit-power over all the 
wireless links. The related analytical framework is available, 
by setting ft = 1, in © and (OS for SSK and TOSD-SSK 
modulations, respectively; ii) r pm = 1 for all the analyzed 
scenarios; iii) the bpcu of SSK and TOSD-SSK modulations 
are equal to R = log 2 (Nt); iv) as far single-antenna and 
Alamouti schemes are concerned, we consider QAM with 
constellation size M and bpcu equal to R = log 2 (M); v) as 
mentioned in Section [V] the shaping filters are obtained from 
[25]. For example, when Nt = 4, w tl (•) in AppendixQ]is used 
for SSK modulation, single-antenna, and Alamouti schemes, 
while the set of four orthogonal filters in (|35T > is used for 
TOSD-SSK modulation. Furthermore, for a fair comparison 
among the modulation schemes, the same spectral efficiency 
(measured in bpcu) is considered; vi) the ABEP of the P- 
CSI scenario is computed by assuming an infinite number of 
pilot pulses; vii) the ABEP of single-antenna and Alamouti 
schemes is obtained through Monte Carlo simulations only, 
but to check that our simulator is well-tuned the numerical 
results are compared, for the P-CSI scenario, to the ABEP 
predicted by the union-bound recently developed in [50] for 
single- and multi-user systems; viii) as far as single-antenna 
schemes are concerned, E m is the average energy transmitted 
for each information symbol; and ix) as far as the Alamouti 
scheme is concerned, E m is the average energy transmitted 
for each information symbol from the two active transmit- 
antennas, i.e., E m is equally split between the two antennas. 
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E m /N [dB] 



Fig. 2. ABEP of SSK modulation against E m /N for: i) N t = 2 (1 bpcu); 
ii) N r = {1, 2,4}; iii) N p = {1, 3, 10}; and iv) P-CSI denotes the ABEP 
with no channel estimation errors. Solid lines show the analytical model and 
markers show Monte Carlo simulations. 
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N =16(4 bpcu) 




E m /N [dB] 



Fig. 4. ABEP of SSK modulation against E m /N for: i) N t = 16 (4 bpcu); 
ii) N r = {1,2,4}; iii) N p = {1,3, 10}; and iv) P-CSI denotes the ABEP 
with no channel estimation errors. Solid lines show the analytical model and 
markers show Monte Carlo simulations. 
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Fig. 3. ABEP of SSK modulation against E m /N for: i) Nt = 8 (3 bpcu); 
ii) N r = {1, 2,4}; iii) N v = {1, 3, 10}; and iv) P-CSI denotes the ABEP 
with no channel estimation errors. Solid lines show the analytical model and 
markers show Monte Carlo simulations. 

The results are shown in Figs. [2^6] for transmission tech- 
nologies with no transmit-diversity gain (SSK and QAM), and 
in Figs. rTT-fTOl for transmission technologies with transmit- 
diversity gain (TOSD-SSK and Alamouti). As far as SSK 
and TOSD-SSK modulations are concerned, we observe 
that: i) our analytical frameworks are very accurate and 
asymptotically-tight for all the analyzed scenarios. In par- 
ticular, as expected, they are exact for Nt = 2; ii) there 
is no loss of the diversity order in the presence of channel 
estimation errors. Only a loss of the coding gain can be 
observed for all MIMO setups; iii) even though in space 
modulation the information is encoded into the CIRs, the 
performance degradation observed when reducing the number 
of pilot pulses, N p , is not very high, and the ABEP is very 
close to the P-CSI lower-bound, in the analyzed scenarios, 
for N p = 10; iv) SSK and TOSD-SSK modulations have 



a SNR penalty, with respect to the P-CSI lower-bound, of 
approximately 3dB and 2dB when N p = 1, respectively; v) 
the ABEP gets worse for increasing Nt, as a consequence 
of the increased size of the spatial-constellation diagram, and 
gets better for increasing N r , due to the receive-diversity gain; 
and vi) TOSD-SSK modulation significantly outperforms SSK 
modulation, due to the transmit-diversity gain introduced by 
the orthogonal pulse shaping design. 

As far as the performance comparison with single-antenna 
and Alamouti schemes is concerned, the following conclu- 
sions can be drawn (see Table HI] for numerical values): i) 
SSK modulation outperforms single-antenna QAM, in all the 
analyzed scenarios, for spectral efficiencies greater than 2 bpcu 
and for N r > 1. If N r = 1, QAM always outperforms 
SSK modulation; ii) SSK and single-antenna QAM have 
almost the same robustness to channel estimation errors, with 
a SNR penalty, with respect to the P-CSI lower-bound, of 
approximately 3dB when N p = T, iii) TOSD-SSK modulation 
outperforms the Alamouti scheme with QAM, in all the ana- 
lyzed scenarios, for spectral efficiencies greater than 2 bpcu. 
In particular, unlike SSK modulation, TOSD-SSK modulation 
is superior to the Alamouti scheme with QAM for N r = 1 as 
well. This is due to the transmit-diversity gain of TOSD- 
SSK modulation; iv) TOSD-SSK modulation is more robust 
to channel estimation errors than the Alamouti scheme with 
QAM. A clear example can be observed in Table [II] when 
R = 2 bpcu and N p = 1. In fact, the Alamouti scheme is 
superior to TOSD-SSK modulation in the P-CSI scenario, 
but TOSD-SSK modulation provides better performance if 
N p = 1 and N r > 1. More in general, Table [IT] shows 
that the Alamouti scheme with QAM has a SNR penalty, 
with respect to the P-CSI lower-bound, of approximately 3dB 
when N p = 1, while TOSD-SSK modulation has a SNR 
penalty of only 2dB; and v) the performance gain of SSK and 
TOSD-SSK modulations with respect to single-antenna and 
Alamouti schemes increases with N r , because, as analytically 
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M = 8 (3 bpcu) 

10° F ! ! ! ! 




E m /N [dB] 

Fig. 5. ABEP of QAM against E m /N for: i) M = 8 (3 bpcu); ii) 
N r = {1,2,4}; iii) N p = {1,3,10}; and iv) P-CSI denotes the ABEP 
with no channel estimation errors. Solid lines with markers or just markers 
show Monte Carlo simulations. Dashed lines show the union-bound computed 
from [50] with no channel estimation errors at the receiver (P-CSI scenario). 
This union-bound is shown only for a subset of curves in order to improve 
the readability of the figure, and avoid overlap among closely-spaced curves. 



than 2 bpcu. In all the cases, the price to be paid for this 
performance improvement is the need of increasing the number 
of radiating elements N t at the transmitter, while still retaining 
a single-RF chain and avoiding inter-antenna synchroniza- 
tion, which are beneficial for low-complexity implementations 
[42]. This remark is somehow similar to [51], as far as the 
achievable transmit-diversity of STBCs is concerned. Finally, 
it is worth emphasizing that the need of a large number 
of radiating elements seems not to be a critical bottleneck 
for the development of the next generation cellular systems, 
as current research is moving towards the utilization of the 
millimeter-wave frequency spectrum [52]. In fact, in this band 
compact horn antenna-arrays with 48 elements and compact 
patch antenna-arrays with more than 4 elements at the base 
station and at the mobile terminal, respectively, are currently 
being developed to support multi-gigabit transmission rates 
[53]. Furthermore, SSK and TOSD-SSK seem to be well- 
suited low-complexity modulation schemes for the recently 
proposed "massive MIMO" paradigm [54], according to which 
unprecedent spectral efficiencies can be achieved in cellular 
networks by using antenna-arrays with very large (with tens 
or hundreds) active radiating elements. 



M = 16(4 bpcu) 




E m /N [dB] 



Fig. 6. ABEP of QAM against E m /N for: i) M = 16 (4 bpcu); ii) 
N r = {1, 2, 4}; iii) N p = {1, 3, 10}; and iv) P-CSI denotes the ABEP with 
no channel estimation errors. Solid lines with markers or just markers show 
Monte Carlo simulations. Dashed lines show the union-bound computed from 
[50] with no channel estimation errors at the receiver (P-CSI scenario). This 
union-bound is shown only for a subset of curves in order to improve the 
readability of the figure, and avoid overlap among closely-spaced curves. 

proved in [50], space modulation takes much better advantage 
of receive-diversity. The results shown in this section confirm 
that this trend is retained in the presence of channel estimation 
errors as well. 

In conclusion, SSK modulation is as robust as single- 
antenna systems to imperfect channel knowledge, and it pro- 
vides better performance when the target spectral efficiency 
is greater than 2 bpcu and N r > 1. On the other hand, 
TOSD-SSK modulation is more robust than the Alamouti 
scheme to imperfect channel knowledge, and it provides better 
performance when the target spectral efficiency is greater 



VII. Conclusion 

In this paper, we have analyzed the performance of space 
modulation when CSI is not perfectly known at the receiver. 
A very accurate and general analytical framework has been 
proposed, and it has been shown that, unlike common be- 
lief, SSK modulation has the same robustness to channel 
estimation errors as conventional modulation schemes, while 
TOSD-SSK modulation is less sensitive to channel estima- 
tion errors than conventional modulations. Also, it has been 
shown that few pilot pulses are needed to achieve almost the 
same performance as the P-CSI lower-bound, and that the 
performance gain, over state-of-the-art MIMO technologies, 
promised by space modulation is retained even with imperfect 
channel knowledge. These results confirm the usefulness of 
space modulation in practical operating conditions, and, in 
particular, the notable performance advantage of TOSD-SSK 
modulation, which provides transmit-diversity and is more 
robust to channel estimation errors than conventional schemes, 
such as the Alamouti code. 
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TABLE II 

Required E m /No (dB) to get ABEP = for all scenarios except SSK modulation and single- antenna QAM if N r = 1, for 

WHICH THE E m /N (dB) TO GET ABEP = 10~ 2 IS SHOWN. FOR SSK MODULATION AND SINGLE-ANTENNA QAM, EACH ROW SHOWS THE E m /N 
(dB) FOR N r = I I N r = 2 I N r = A. FOR SSK MODULATION AND THE ALAMOUTI SCHEME WITH QAM, EACH ROW SHOWS THE E m /N (dB) FOR 

N r = 1 / N r = 2. The values have an error approximately equal to ±0.1dB. 



SSK 


Rate 


N p = l 


N p = 3 


N p = 10 


P-CSI 


1 bpcu 


22.9 1 25.3 / 16.2 


21.1 /23.5 / 14.5 


20.3 / 22.7 / 13.6 


19.9/22.3 / 13.2 


2 bpcu 


26 / 26.8 / 17 


24.2 / 25.1 / 15.4 


23.4 / 24.3 / 14.5 


23 / 23.8 / 14 


3 bpcu 


29 / 28.4 / 17.9 


27.3 / 26.6 / 16.2 


26.4/ 25.8 / 15.3 


26 / 25.4 / 14.9 


4 bpcu 


32 / 29.9 / 18.7 


30.3 / 28.1 / 17 


29.5 / 27.3 / 16.2 


29 / 26.9 / 15.7 


Single-Antenna QAM 


Rate 


N p = l 


N p = 3 


N p = 10 


P-CSI 


1 bpcu 


19.8 / 22.3 / 13.1 


18.2/20.6 / 11.5 


17.2 / 19.7 / 10.6 


16.8 / 19.3 / 10.1 


2 bpcu 


22.7 / 25.2 / 16.2 


21.1 /23.5 / 14.5 


20.3 / 22.7 / 13.6 


19.9/22.2 / 13.2 


3 bpcu 


27.5 / 29.9 / 21.1 


25.6 / 28 / 19.1 


24.9 / 27.4 / 18.2 


24.6 / 27 / 18.1 


4 bpcu 


29.6 / 32 / 23.3 


27.8 / 30.1 / 21.3 


27.1 / 29.6 / 20.5 


26.8 / 29.1 / 20.3 


TOSD-SSK 


Rate 


N p = l 


N p = 3 


N p = 10 


P-CSI 


1 bpcu 


27.2 / 18.2 


26 / 16.9 


25.5 / 16.4 


25.3 / 16.2 


2 bpcu 


28.7 / 19 


27.5 / 17.8 


27 / 17.3 


26.8 / 17 


3 bpcu 


30.2 / 19.8 


29 / 18.6 


28.5 / 18.2 


28.4/ 17.8 


4 bpcu 


31.9 / 20.7 


30.5 / 19.4 


30.1 / 18.9 


29.9 / 18.7 


Alamouti QAM 


Rate 


N p = l 


N p = 3 


N p = 10 


P-CSI 


1 bpcu 


25.3 / 16.2 


23.5 / 14.5 


22.8 / 13.5 


22.3 / 13.2 


2 bpcu 


28.4 / 19.3 


26.5 / 17.5 


25.7 / 16.6 


25.4 / 16.3 


3 bpcu 


32.9 / 24 


31.4/22.2 


30.4 / 21.3 


30 / 21 


4 bpcu 


35.2 / 26.2 


33.3 / 24.3 


32.6 / 23.5 


32.3 / 23.3 



' (0 - (-Ties) * «) + * (0 + G/S) ** (0 + {-^) P, (0 + (-^jj) p 5 (0 

( «*> (0 = (-Tin) * «) + ("^) ft (0 + G/i) *> (0 + {^W) v, (0 + (-^5) p 6 (0 

«*. (0 = G/l) * (0 + (-^P) P2 (0 + (0) P3 (0 + (^^) (0 + (-^) (0 

^ (0 = G/i) w (0 + (^^) p 2 (0 + (o)i* (0 + (-^^) ^4 (0 + (-^) P5 (0 



' wii M = (-Tin) Pl M + f^?^) P2 (") + G/S) P3 M + (-^^) * (-) + (-Trio) p * M 



V^ t3 (w) 



Pi («) + p 2 («) + (0)P 3 M + (^5^) p 4 M + (-^) p 5 («) 

wi. («) = G/&) Pi M + (^±^) p 2 ( W ) + (o) p 3 ( W ) + (-^^j p 4 M + (-^) p 5 («) 



(36) 
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N = 2 (1 bpcu) 




E m /N [dB] 



Fig. 7. ABEP of TOSD-SSK modulation against E m /N for: i) N t = 2 
(1 bpcu); ii) N r = {1,2}; iii) N p = {1,3,10}; and iv) P-CSI denotes 
the ABEP with no channel estimation errors. Solid lines show the analytical 
model and markers show Monte Carlo simulations. 
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N = 8 (3 bpcu) 




E m /N [dB] 



Fig. 9. ABEP of TOSD-SSK modulation against E m /N for: i) N t = 8 
(3 bpcu); ii) N r = {1,2}; iii) N p = {1,3,10}; and iv) P-CSI denotes 
the ABEP with no channel estimation errors. Solid lines show the analytical 
model and markers show Monte Carlo simulations. 



N ( =16(4bpcu) M = 16 (4 bpcu) 




E m /N [dB] E m ' N oI dB ] 



Fig. 8. ABEP of TOSD-SSK modulation against E m /N for: i) N t = 16 
(4 bpcu); ii) N r = {1,2}; iii) N p = {1,3,10}; and iv) P-CSI denotes 
the ABEP with no channel estimation errors. Solid lines show the analytical 
model and markers show Monte Carlo simulations. 



Appendix I 
Orthogonal Shaping Filters for N t = 4 



Fig. 10. ABEP of Alamouti scheme with QAM against E m /No for: i) 
M = 16 (4 bpcu); ii) N r = {1,2}; iii) N p = {1,3,10}; and iv) P-CSI 
denotes the ABEP with no channel estimation errors. Solid lines with markers 
or just markers show Monte Carlo simulations. Dashed lines show the union- 
bound computed from [50] with no channel estimation errors at the receiver 
(P-CSI scenario). This union-bound is shown only for a subset of curves 
in order to improve the readability of the figure, and avoid overlap among 
closely-spaced curves. 



In this appendix, we show an example of orthogonal shaping 
filters that can be used for TOSD-SSK modulation. Without 
loss of generality we consider the case study with Nt = 4, but 
the procedure can be generalized to larger antenna-arrays. 

More specifically, we consider the procedure described 
in [25], which allows us to generate orthogonal shaping 
filters with the same time-duration and bandwidth. Similar 
techniques are available in [46], [47]. From [25], we can 
obtain the four orthogonal impulse (time) responses shown 
in d35l > at the bottom of the previous page, as well as the 
four related frequency responses (Fourier transform) P (w) = 
(1 /V2n) / +~ p (£) exp (-jug) d£ shown in (f36]l at the bot- 



tom of the previous page too, where we have defined: 
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and: 



Pi (w) 




12to-48*oW 2 + 16tow' 1 



■5 (*<><*>) 



exp 



I (*ow) 5 



(38) 

Finally, we mention that, by adjusting the form factor to, 
the bandwidth can be arbitrarily chosen, and both narrow- and 
wide-band communication systems can be considered. 
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